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STM observation of electronic wave interference effect in finite-sized graphite 
with dislocation-network structures 
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Superperiodic patterns near a step edge were observed by Scanning Tunneling Microscopy on several-layer- 
thick graphite sheets on a highly oriented pyrolitic graphite substrate, where a dislocation network is generated 
at the interface between the graphite overlayer and the substrate. Triangular- and rhombic-shaped periodic pat- 
terns whose periodicities are around 100 nm were observed on the upper terrace near the step edge. In contrast, 
only outlines of the patterns similar to those on the upper terrace were observed on the lower terrace. On the 
upper terrace, their geometrical patterns gradually disappeared and became similar to those on the lower terrace 
without any changes of their periodicity in increasing a bias voltage. By assuming a periodic scattering poten- 
tial at the interface due to dislocations, the varying corrugation amplitudes of the patterns can be understood as 
changes in the local density of states as a result of the beat of perturbed and unperturbed waves, i.e., the inter- 
ference in an overlayer. The observed changes in the image depending on an overlayer height and a bias voltage 
can be explained by the electronic wave interference in the ultra thin overlayer distorted under the influence of 
dislocation-network structures. 

PACS numbers: 68.37.-d, 68.37.Ef, 72.10.Fk,73.90.+f 



I. INTRODUCTION 

Scanning tunneling microscopy (STM) observations of su- 
perperiodic patterns on metal surfaces have been reported in 
several finite-sized systems. They are ascribed to interference 
patterns of free electron waves scattered by adatoms and step 
edges, for example, an Ag(lll) surface near a step edge, 7 777 
a Cu(l 1 1) surface surrounded by 76 Fe adatoms, ? 7 7 and so 
on. These reports have clarified that scattered and interfered 
waves on the surface can be observed as periodic patterns 
which are related to the Fermi surface of bulk and surface 
states of metals where free electrons can move around. Re- 
cently, a superperiodic pattern has been also reported in semi- 
conductor surfaces such as InAs/GaAs(l 1 1)A. 7 7 In this case, 
the pattern is also an interference pattern, which is gener- 
ated by electron waves scattered at step edges on semicon- 
ductor surfaces because a two-dimensional (2D) electron gas 
is generated due to the band bending by the surface recon- 
struction. This phenomenon is interesting and characteristic 
of the surface electronic structure of isotropic semiconduc- 
tors; that is, generated electrons whose characters resemble 
free electrons in metals can move in a surface thin layer in 
spite of three dimensionality in the electronic structures of 
semiconductors. As for 2D electronic systems, the present 
authors have observed the electronic wave interference effect 
on nanographene sheet inclined with respect to a highly ori- 
ented pyrolytic graphite (HOPG) substrate by STM. ? Here, 
a nanographene sheet interacts very weakly with the HOPG 
substrate, where electrons are confined in the 2D sheet and 



the in-plane potential changes gradually. 

Meanwhile, superperiodic patterns on an HOPG sub- 
strate observed by STM have been also reported in many 
papers. 7 ???????? Those patterns are not generated by elec- 
tronic wave interference effects, but are caused by multiple- 
tip effects, rotational stacking faults, and dislocation-network 
structures. A pattern caused by multiple-tip effects originates 
from superimposing two different information of graphite lat- 
tice in one domain imaged by a tip apex and that in another 
domain, with a relative rotation, imaged by a metal contami- 
nation that is attached to the tip. A moire pattern and a pattern 
caused by the dislocation-network structures result from the 
spatially varied local density of states (LDOS), which are re- 
lated to the stacking faults by the relative rotation between 
two adjacent graphene layers and by the lattice distortion at 
the interface, respectively. Graphite with a stacking fault can 
be represented as abcab . . . where represented as ababa . . . 
for an ordinary stacking of graphite and c for a faulted layer. 
The periodicity of a moire pattern and a pattern caused by 
the dislocations can be explained by an angle of the relative 
rotation 7 7 7 7 and by the periodic domain of stacking faults 
generated in a slip plane, respectively. Among the reports, a 
change in the bias voltage interestingly induces a change in 
the periodicity of the superperiodic patterns that come from 
the dislocation network, similar to that is observed by trans- 
mission electron microscopy (TEM). 7 

LDOS calculation of faulted stacking parts cannot repro- 
duce the corrugation amplitude of superperiodic structures 
that have been reported so far. 7 7 7 About a moire pattern and 
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a dislocation-induced pattern, one can find that the interface 
between an overlayer and a substrate is taken as a scattering 
layer and that the overlayer is regarded as a finite-sized re- 
gion in the normal direction to the surface. 7 Electron waves 
normal to the surface can be scattered by the surface and the 
interface, resulting in the generation of standing waves. This 
is the electron confinement effect in the ID direction normal 
to the surface, which is very important for the superperiodic 
LDOS at the surface in terms of the corrugation amplitude in 
a STM image. The corrugation amplitude of a superperiodic 
pattern is expected to depend on a bias voltage of STM and the 
overlayer thickness as the character of waves in the overlayer. 
In this paper, we report on the observation of different super- 
periodic patterns that originate from the dislocation-network 
structures, on both terraces near a step edge and present their 
bias voltage dependence of the corrugation amplitudes with 
no change in the periodicity. The patterns are explained as 
the spatially varied LDOS affected by the interference in the 
overlayer. 



II. EXPERIMENTAL 

All images in the present paper were observed by using 
a commercial STM system (Digital Instruments, Nanoscope 
E) under an ambient condition at room temperature with the 
constant-current mode at 0.7 nA using a mechanically cut Pt- 
Ir tip. Sample bias voltages for these observations were var- 
ied from near the Fermi level to higher voltages, typically at 
0.02, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, and 0.6 V, except for that in 
Fig. 2(b). The sample was fabricated by the heat-treatment of 
an HOPG substrate at 1600 °C in Ar flow after cleaving it by 
an adhesive tape for obtaining a fresh surface. It is possible 
that dislocations were generated at several layers beneath the 
surface during the heat-treatment process. 



III. RESULTS 

Superperiodic patterns were observed on the surface by 
STM, which extended over several /im 2 , and a part of the 
area is shown in an 1.5x1.5- fita 2 image in Fig. 1(a), where 
a bias voltage was 0.2 V. In this image, there are triangular- 
, rhombic-, and net-shaped patterns whose periodicities are 
around 100 nm but gradually changed depending on the posi- 
tion. There are also complicated patterns that seem to be the 
superimposed of those patterns. The lines pointed by arrows 
a denote step edges of graphite, while the line pointed by ar- 
rows (3 denotes a domain boundary where the difference in 
the heights between the two regions faced at the boundary is 
much less than the interlayer distance of graphite (0.335 nm 
in the bulk). The presence of the step edge and the domain 
boundary is confirmed by the cross-sectional profiles shown 
in Fig. 1(a). The lower terrace in Fig. 1(a) extends to the left 
direction by about 3.5 //m and is terminated by a boundary 
between the graphene overlayer and the HOPG substrate. Pat- 
terns of those shapes change into y-shaped and linear patterns 
near the graphene overlayer edge and end at the edge (the re- 




FIG. 1: (a) STM image (1.5x1.5 //m 2 ) of superperiodic patterns ob- 
served at V,=0.2 V. There are triangular-, rhombic- and net-shaped 
patterns whose periodicity is around 100 nm. There are also com- 
plicated patterns in some parts, where two types of patterns are su- 
perimposed. Arrows a and /? denote the positions of a step edge of 
graphite and a domain boundary, respectively. Cross-sectional pro- 
files are taken along the horizontal lines, P-l to P-3, whose lengths 
are 500 nm. (b) Magnified STM image (500x500 nm 2 ) of the center 
region of (a), showing two types of patterns. The heights of the lower 
and upper terraces correspond to two and three graphene layers from 
the substrate, respectively. 
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gion showing the y-shaped and linear patterns is not shown). 
The triangular- and net-shaped patterns similar to those in 
Fig. 1 (a) have been previously observed in STM, ? 7 TEM, ? 7 7 
and other investigations. 7 7 7 In those reports, diffraction pat- 
terns in TEM images and superperiodic patterns in STM im- 
ages were attributed to the modified LDOS caused by rhom- 
bohedral stacking faults due to partial dislocations. The partial 
dislocations are defined by the Burgers vector that converts an 
afe-stacked layer in ordinary graphite to an ac-stacked layer 
with respect to a glide plane. The conversion of stacking oc- 
curs abruptly accompanied with a lattice distortion where a 
sharp-edged periodic pattern is generated. Therefore, the pat- 
terns in Fig. 1(a) are considered to come from dislocations at 
the interface between the graphite overlayer and the HOPG 
substrate from the shapes and the average periodicity. 

A magnified image of Fig. 1(a) near a step edge is shown in 
Fig. 1(b). The height of the lower terrace at the bottom left in 
Fig. 1(b) from the HOPG substrate is 0.67+0.02 nm in average 
(for V s =0.2 V) from the cross-sectional profile analyses of the 
observed image at the boundary between the graphene over- 
layer and the substrate (not shown). The value corresponds to 
a thickness of two graphene layers from the substrate. The 
upper terrace at the center part in Fig. 1(b) has a height of 
three graphene layers from the substrate, as estimated from 
the cross-sectional profile of the step edge, whose height dif- 
ference is 0.39-0.41 nm (for V s =0.2 and 0.5 V) including 
the corrugation amplitude of superperiodic patterns [P-l in 
Fig. 1(a)]. The image at a low bias voltage of 0.02 V near the 
Fermi energy is shown in Fig. 2(a). Though this image was 
obtained at almost the same place as that shown in Fig. 1(b), 
there are a few differences in the contrast and the shape of the 
patterns. In Fig. 2(a), three regions are indicated; regions A, B, 
and C contain a triangular-shaped pattern, a rhombic-shaped 
pattern, and a net-shaped pattern, respectively. In intermedi- 
ate regions A-B and B-C, there are complicated contrasts that 
are superimposed of patterns in two regions. In regions A 
and B, the apparent height of lines which divide the patterns 
into individual geometric units was lower than the center of 
the unit by about 0. 1 nm. Crossed points of lines were fur- 
ther depressed from the lines by about 0.1 nm, resulting in 
the "contracted nodes" in the image. In region C, however, 
lines are imaged higher than the center of the unit by about 
0.05 nm and crossed points of the lines are the highest (about 
0.005 nm higher than the lines), giving the "extended nodes" 
in the image. Except for the slight contrast, the patterns in 
regions B and C appear to have contrast inverted from each 
other. In Fig. 1(b), however, the part corresponding to region 
B is not the same pattern as region B in Fig. 2(a), suggesting 
that the corrugation amplitudes and the shapes of patterns de- 
pend on a bias voltage. The triangular shape in region A is 
almost the same as can be seen in the comparison of Fig. 1(b) 
and Fig. 2(a). However, the contracted nodes on the upper ter- 
race at a bias voltage of 0.02 V [Fig. 2(a)] changed to the ex- 
tended nodes at 0.2 V [Fig. 1(b)]. Figure 2(b) is a magnified 
image near a contracted node of the upper terrace in Fig. 2(a), 
which is marked by a black dot, taken at V s =0.002 V, 7=1.7 
nA. A straight line drawn on triangular lattice points at the 
bottom right part is extended to the valley sites of the trian- 




FIG. 2: (a) STM image (500x500 nm 2 ) of the superperiodic patterns 
at a low sample bias voltage of 0.02 V, which is expected to reflect 
the density of states close to the Fermi level. The imaged area is 
shifted from that of Fig. 1 to the bottom by about 200 nm. A, B, and 
C denote regions of triangular-, rhombic- and net-shaped patterns, 
respectively. Arrows indicate complicated patterns, where two pat- 
terns are superimposed. Lines that divide the geometric patterns into 
individual units cross at contracted nodes in regions A and B, and 
at extended nodes in region C. The apparently depressed contrast 
neighbors to the lines are artificial effect to make the image clearer 
in region C. (b) Atomically resolved STM image (6.0x6.0 nm 2 ) of 
one individual triangular pattern near a contracted node on the upper 
terrace in (a), which is marked by a black dot, at Kt=0.002 V and 
7=1.7 nA. A straight line placed on triangular lattice sites at the bot- 
tom right part is extended to the valley sites of the triangular lattice 
at the top left part through a distorted lattice part. 
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(a) (b) 




FIG. 3: STM images (500x500 nm 2 ) of superperiodic patterns at higher sample bias voltages; (a) V s =03 V, (b) V,=0.4 V, and (c) V s =0.5 V. 
By increasing the bias voltage, the corrugation amplitude of superperiodic patterns on the upper terrace decreased gradually [(a) and (b)] and 
changed into a net pattern (c). In contrast, no significant change was observed for the pattern on the lower terrace. The net pattern appearing 
on the upper terrace of (c) is similar to that on the lower terrace. Height differences between two points depicted in (a) are shown in (d) for 
clarifying the bias-dependent contrast. Solid and blank circles are the height differences of the upper and lower terrace, respectively. (Circles 
at the sample bias of around V are the height differences at V,=0.02 V.) 



gular lattice at the top left part, indicating the presence of a 
distortion at the center part of the image. This atomically re- 
solved image supports that the observed patterns come from 
the dislocation-network structure. As for patterns at higher 
bias voltages, Figs.3(a)-(c) display images of almost the same 
places as Figs. 1(b) and 2(a) at bias voltages of 0.3 V, 0.4 V, 
and 0.5 V, respectively. The patterns in Fig. 3(a) seem to re- 
semble those in Fig. 1(b), except for the pattern change from 
rhombic-shaped to triangular-shaped in region B. In regions 



A and C, the corrugation amplitudes of patterns on the upper 
terrace in Fig. 3(a) are smaller by about 0.03 nm than those 
in Fig. 1(b), whereas those on the lower terrace in Fig. 3(a) 
are larger by about 0.06 nm than those in Fig. 1(b), as shown 
in Fig. 3(d). For clarity, the height differences between two 
points on each terrace dependent on the bias voltage is shown 
in Fig. 3(d). Just by increasing a bias voltage, patterns on the 
upper terrace in Fig.3(a) are changed into a net-shaped pat- 
tern on the lower terrace in Fig. 3(c). They are similar to the 
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pattern on the lower terrace in Figs.3(a)-(c), however, slightly 
inversed contrast is observed as shown in Fig. 3(a) similar to 
that in region C in Fig. 2(a). Changes of corrugation ampli- 
tudes, with a maximum at V. s =0.3 V are observed in the pat- 
terns on the lower terrace as shown in Fig. 3(d). A similar 
image to Fig. 3(c) was also observed at V s =0.6 V (not shown). 
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IV. THEORETICAL MODEL AND DISCUSSION 

According to previous reports, the periodicity of super- 
periodic patterns changed dependent on a bias voltage or 
just by scanning the tip, which was attributed to dislocation 
motion. 7 ? However, the periodicity of the observed patterns 
in the present study did not change in the range of voltages 
used for imaging (0.02 - 0.6 V). Therefore, the observed 
phenomenon is different from that in Ref. 15 and 16, which 
show a dislocation motion induced by the applied bias volt- 
age. A dislocation motion is not generated in the observed 
phenomenon in the present paper. Instead, only the superpe- 
riodic corrugation amplitudes of the observed patterns varied, 
without any change in the periodicity, depending on an over- 
layer height from the substrate and a bias voltage of STM, 
as shown in Fig. 3(d). This is the first observation of the 
bias-dependent contrast and pattern shapes of superperiodic 
patterns on graphite without any change of the periodicity.// 

First, we discuss the origin of the bias-dependent corruga- 
tion amplitudes. As shown in Figs. 1-3, it appears that the 
patterns on the upper and lower terraces are connected con- 
tinuously at the step edge, independent of the bias voltage. 
This suggests that the patterns observed on both terraces come 
from the same origin. It should be noted that the experimen- 
tal results cannot be explained simply by calculating the DOS 
of faulted stacking, because the observed patterns at differ- 
ent terraces have contrast inverted from each other as shown 
in Figs. 1 -3(b). It seems natural to assume that an array of 
faulted stacking is not changed abruptly across the step edge 
if the dislocation network is continuous at the interface. We 
cannot also explain the property, on the basis of the faulted 
stacking, that the superperiodic corrugation amplitudes on the 
lower terrace become larger although the gap between the tip 
and the sample becomes larger in increasing a bias voltage 
from 0.02 to 0.3 V as shown in Fig. 3(d). In other words, the 
observed behavior is considered to be due to the LDOS at the 
surface, taking into account the fact that observed corruga- 
tion amplitudes on the lower terrace become larger. Then, the 
LDOS should explain the gradual decreases of the corrugation 
amplitudes and the variations of patterns on the upper terrace 
in increasing the bias voltage without changing the periodic- 
ity of the patterns, and that should also explain the increase of 
the corrugation amplitudes on the lower terrace in increasing 
a bias voltage near the Fermi energy. 

Here, we discuss the interference effect for explaining the 
bias-voltage dependence of superperiodic patterns on the ba- 
sis of a theoretical treatment reported in Ref.21. Consider- 
ing the scattering potential at the interface, one can find that 
the LDOS at the surface is related to the interference effect 
of electrons that are scattered at the surface and the interface 
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FIG. 4: Model of the potential at the interface between the graphite 
overlayer and the substrate, (a) The cross-sectional profile of square 
potential along the x or y axis. The periodicity is 2L and the potential 
is L/3 in width and 2voS(z) in height, (b) The projection of the square- 
patterned potential on the xy plane. Gray lines represent potential 
lines and black squares represent potential nodes, whose potential 
height is the sum of the ID potentials in x and y axes. A height of the 
potential nodes is twice as large as that of the potential lines, (c) The 
position of the surface and the interface along the z axis. The surface 
and the interface are located at I and in the z axis, respectively. 



between the overlayer and the substrate. The LDOS at the 
surface can be given as sin 2 (fe) using coordinate z and a wave 
number k along the axis normal to the surface (the z axis) in 
case that the lateral wave number of a superperiodic pattern 
nearly equals to by comparison with the wave number orig- 
inating from the lattice. If one treats a scattering potential 
at the interface by perturbation, a beat can be generated by 
the interference between the perturbed and the unperturbed 
waves. In this case, the LDOS at the surface is proportional to 
sin(fe) cos(fe'z), where k' and k are a perturbed wave number 
and an unperturbed wave number, respectively. Next we will 
show the detailed derivation. 

In the present discussion, a square-patterned potential with 
a periodicity of 2L at the interface is employed as shown in 
Fig. 4, for a calculation of the probability density of the wave 
function confined in the plane for generating an abrupt poten- 
tial change associated with the dislocation-network structures. 
We place a square potential with L/3 in width and 2vq6(z) in 
height, where L is the half of the periodicity of the square po- 
tential and vq is the strength of the scattering potential, at the 



6 



line dividing the patterns into geometrical units as a simple 
model to reproduce the patterns in regions B and C, as shown 
in Fig. 4(a) and (b). Though the square-shaped pattern in the 
present model is different from the experimental result (the 
rhombic-shaped pattern in regions B and C, the triangular- 
shaped pattern in region A, and the complicated pattern in 
their intermediate regions), it can make a theoretical treatment 
easier with any loss of validity. (Note that the problem of the 
square-patterned potential in rectangular coordinate can be re- 
duced to the ID problems along the x and y axes.) The repro- 
duction of a pattern in region A is beyond the present model 
since the shape of a pattern in region A is complicated to solve 
in the similar manner that is applied for regions B and C. If we 
locate the surface and the interface positions at I and 0, respec- 
tively, in the z axis as shown in Fig. 4(c) and introduce the delta 
function 6(z) at the interface, this potential can be expressed 
using the Fourier analysis 

n 

x (e i<? '" * + e _i<? " ' x + e 1 *" + e _i *" y ), ( 1 ) 

where ti is the Planck constant over 2n, m ± is the effective 
mass along the z axis, an is the nth component which equals 
to {2/(n7r)}{sin(«7r) - sin(5«7r/6)} for the square potential, and 
q xn and q yn , which take discrete values (nn/L) (n-1,2, . . .), 
are the nth wave vectors in the x and the y axes, respectively. 
Assuming the linear combination of in-plane plane waves and 
wave function A q ^{z) for the z component, the wave function 
is represented to be 

^W)^^*^''' 1 (2) 

?*.?>. 

Based on this wave function and the connecting condition de- 
rived from the Schrodinger equation with the square-patterned 
potential, the perturbed wave function is given as 

^(Jcy.z) = v Q J]{a n /(ik')}(l - e 2lW ) 

n 

X (e'*' Z - e - i *'z+ 2ii '') e ± (?-m ,JC +?ra-)') (3) 

where the term of v 2 , is neglected because of its small contri- 
bution. For the unperturbed wave, we take a plane wave in the 
direction of the z axis, 

= e lfe , (4) 

with 

k' 2 = k 2 -(m n /m ± )\q n \ 2 , (5) 

where my is the effective mass in the xy plane. In the over- 
layer shown in Fig.4(c), the unperturbed waves can coex- 
ist with the perturbed waves. Since the total wave function 
^totaiC*! y, z) is the sum of y ¥+, *P_, and *Fo, the probability 



density, | ^totaiC^.y.z) I 2 is represented as follows: 

|TW*,y,z)| 2 =4sin 2 (fe) 

- 32 y^(a n vo/fc / ) sin(£/) cos(£7) 

n 

X sin(fe) sin(fc'z) 

x {cos(q m • x) + cos{q yn • y)}. (6) 

Here, the xy plane is shifted from (z - 1) to z for simplicity and 
the very small term of v 2 , is also neglected. Using the Tamm 
states near the surface (z ~ 0), k cot(fcz), and k' cot(fe'z) can be 
regarded as the constant numbers, provided that kz ~ and 
k'z ~ 0. Therefore, the probability density around the surface, 
I ^PtotaiC*. y> 0) | 2 , is roughly expressed with an overlayer height 
from the substrate, I, the perturbed wave number, k', and the 
unperturbed wave number, k: 

I ^totaiOt.y.O) | 2 = -c^^a n ksm(kl)cos(k'l) 

n 

X {cos(q xn ■ x) + cos(<7 > .„ • y)} + const., (7) 

where c and the second term are positive constants and the 
second term is larger absolute value than that of the first 
term. In this equation, the spatially varied probability density, 
which gives a superperiodic pattern, corresponds to the term, 
{cos(q xn -x) + cos(g y „ -y)}. By attributing the unperturbed wave 
to the wave function in bulk graphite, the energy dispersion 
can be given using parameter m ± , the interlayer distance c, 
and the interlayer resonance integral, yi(=0.39 eV) (Ref.28): 

E = ti 2 k 2 /(2m A _)-2y u (8) 

with 

m ± = « 2 /(2c 2 ri ). (9) 

The investigation of the spatially varied LDOS is important in 
order to look over the contrast image of STM from the corru- 
gation amplitude of a superperiodic pattern that depends on a 
bias voltage. In this connection, the difference of the LDOS at 
the surface, {| ^,(0,0,0) | 2 - | ^^(^,0,0) | 2 }, where the 
former and latter terms represent the LDOS at the center and 
edge of an individual geometric pattern unit, respectively, can 
give a simple diagnosis in mapping a superperiodic pattern 
because a potential height is constant except for the edge part 
with a fine oscillation resulted from the Fourier analysis. In 
Fig.5, the difference of the LDOS in an arbitrary unit varies 
as a function of an overlayer height from the interface and 
a bias voltage. It obviously changes in the present range of 
an overlayer height from the interface and a bias voltage that 
were used for STM observation. The positive value of the dif- 
ference means that the LDOS at the bias voltage is larger at 
the center position, (x,y)=(0, 0), surrounded by four potential 
lines than on the potential line, (x,y)=(L, 0), giving a square- 
shaped pattern. Conversely, the negative value suggests that 
the LDOS at the center position is smaller than that on the po- 
tential line, the superperiodic pattern being a net-shaped pat- 
tern. As the absolute value of the difference is proportional 
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Overlayer height from the interface (layer) 

FIG. 5: The difference of LDOS at two points, 

{I 'PtotaiCO.O.O)! 2 - | IVtaifL^O) | 2 }, as a function of an over- 
layer height from the interface and a bias voltage. The periodicity 
of a square potential, 2L, is 70 nm. The difference of the LDOS 
is shown in scale bar and its unit is arbitrary. When the difference 
is a positive (negative) value, a calculated superperiodic pattern 
is a square-shaped (net-shaped) pattern. The unit number of the 
overlayer height (the number of layers) corresponds to the interlayer 
distance between adjacent graphene layers (0.335 nm). 



to the superperiodic corrugation amplitudes, the increase, de- 
crease, and inversion of a corrugation amplitude of the super- 
periodic pattern can be generated depending on a bias voltage 
and an overlayer height from the substrate, as shown in Fig. 5. 
For a comparison between the theoretical model and the ex- 
perimental results, an overlayer height from the interface, /, is 
given in the unit of a single graphene layer thickness, which 
corresponds to the interlayer distance of bulk graphite (0.335 
nm). Assuming that the interface is located at the intermediate 
plane between the overlayer and the substrate (an ideal inter- 
face is shifted by a half of monolayer from the substrate), / is 
given as 



I = s - 0.5 + A, 



(10) 



where s is an overlayer height from the substrate in the same 
unit as that of /, and A is a fitting parameter for a wavy struc- 
ture of graphite. 7 

Figure 6 shows the calculated LDOS in a 2L x 2L square of 
the individual geometrical pattern unit at different overlayer 
heights and varied bias voltages. The extended and contracted 
nodes appear at the crossing points (x,y)=(±L, +L),(+L, +L). 
At 1=1. 5(5=2.0, A=0), the difference of the LDOS enhances 
with increasing a bias voltage from 0.02 to 0.3 V, where the 
calculated superperiodic patterns are net-shaped patterns with 
the extended nodes as clearly seen in Figs.6(a)-(d). This result 
agrees with the experimental evidence that the inverted super- 
periodic patterns (net-shaped patterns) were observed on the 
lower terrace (two graphene layers high from the substrate) 
near the step edge, as shown in Figs. 1-3. Indeed, the increase 
of the corrugation amplitude in Figs. 2(a) and 3(a) can be un- 
derstood on the basis of the increased difference of the LDOS 
because the corrugation amplitude is roughly proportional to 
the LDOS, as mentioned above. At /=2.7(s=3.0, A=0.2), the 
difference of the LDOS decreases in increasing a bias voltage 
from 0.02 to 0.4 V, where the calculated superperiodic pat- 
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FIG. 6: The calculated LDOS in a 2L x 2L square (upper) and its 
cross-sectional profile passing through the center of the geometri- 
cal unit (lower) at different overlayer heights and bias voltages: (a) 
/=1.5 (layer), V,=0.02 (V), (b) 1=1.5, V,=0.30, (c) /=1.5, V s =0.40, 
(d) /= 1.5, V,=0.50, (e) 1=2.1, V s =0.02, (f) 1=2.7, V>0.30, (g) 1=2.7, 
V s =0.40, and (h) /=2.7, V s =0.50. (Top pictures) Lighter brightness 
indicates a higher LDOS value. (Bottom pictures) The x or y value 
in the lateral axis (unit: nm), the LDOS in the vertical axis (arbitrary 
units). 



terns become square-shaped patterns with contracted nodes 
as shown in Figs.6(e)-(g), although the pattern shape (square 
shaped) is different from that of experimental results (rhom- 
bic and triangular shaped). By further increase in the bias 
voltage (0.4 to 0.5 V), the difference of the LDOS has a neg- 
ative value with a net-shaped pattern with extended nodes in 
Fig. 6(h). This result supports the experimental evidence that 
the corrugation amplitude of the superperiodic pattern on the 
higher terrace (three graphene layers high from the substrate) 
decreases gradually and that the superperiodic pattern changes 
into an inverted pattern with increasing a bias voltage further, 
as shown in Fig. 3. Eventually, the observed superperiodic pat- 
terns can be explained by the dislocation network at the in- 
terface and an interference in the overlayer dependent on its 
thickness and bias voltages. 

However, this model cannot explain the relation between 
patterns in regions A and B, and the change in the corrugation 
amplitude of extending and contracting nodes. The potential, 
independent of the bias voltage in the present model, may be 
responsible for the discrepancy. A more appropriate way of 
faulted stacking or a slight relaxation is expected to improve 
the model, including a change of pattern shapes. The simple 
model used in this paper suggests that nodes can be alternated 
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between the extended and the contracted due to the LDOS 
affected by the interference in the overlayer without chang- 
ing the way of stacking. Other remaining problems are the 
discrepancy between the semimetallic electronic structure of 
graphite and the present theoretical model, and a phase shift 
of electron waves due to the Coulomb repulsion in the over- 
layer. These problems call for further investigations in the 
near future. 



and a pattern at the surface can be affected by the beat of then- 
waves. On the basis of the free electron model with the effec- 
tive mass, the corrugation amplitudes of the patterns, which 
are related to the LDOS, are found to vary depending on a 
bias voltage and an overlayer height, and the changed corru- 
gation amplitude of a superperiodic pattern can be understood 
as a change of the LDOS originating from the interference in 
the overlayer. 



V. CONCLUSION 

Superperiodic patterns that come form the dislocation- 
network structure have been observed by STM, where the 
shape and the corrugation amplitudes change dependent on a 
bias voltage and an overlayer height from the substrate with- 
out any variation of their periodicity. Near a step edge, the 
dislocation network that causes patterns on the upper and the 
lower terraces seems to be continuous. By assuming the same 
scattering potential at the interface between both terraces, a 
perturbed wave that generates a superperiodic pattern in the 
plane and an unperturbed wave can interfere in the overlayer, 
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